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一 类 具 时 滞 和 收获 的 捕食 模型 的 稳定 性 与 Hopf 分 支 * 


田 晓 红 ， 徐 瑞 ， 王 丽 丽 
(军械 工程 学 院 应 用 数学 研究 所 ， 石 家 庄 050003) 

摘 要 : 本 文 研 究 一 类 具有 常数 收获 率 和 时 滞 的 捕食 模型 ， 其 中 时 滞 描 述 了 捕食 种 群 的 妊娠 期 。 通 过 分 析 
特征 方程 ， 得 到 了 正平 衡 点 局 部 稳定 的 条 件 。 当 时 滞 了 增加 时 ， 正 平衡 点 失去 稳定 性 ， 当 T 跨 过 
临界 值 时 系统 将 出 现 Hopf 分 支 。 应 用 中 心 流 形 定理 和 规范 型 理论 ， 得 到 了 确定 Hopf 分 支 方向 和 
分 支 周 期 解 的 稳定 性 的 计算 公式 。 最 后 对 所 得 理论 结果 进行 了 数值 模拟 。 
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种 群生 态 学 是 生态 学 中 的 一 个 重要 分 支 ， 也 是 与 人 们 的 生产 生活 最 密 不 可 分 的 学 科 之 一 。 
人 们 对 种 群生 态 学 的 研究 ， 一 方面 要 对 种 群 的 发 展 变化 有 定量 的 分 析 和 预测 ， 另 一 方面 也 是 最 
为 重要 的 就 是 建立 具有 收获 率 的 种 群 系统 ， 通 过 对 系统 的 分 析 ， 判 断 出 付出 多 大 的 捕获 量 ， 才 
可 以 既 维 持 生 态 系统 的 平衡 ， 又 能 使 收获 量 达到 最 大 ， 以 满足 人 类 的 需要 。 因 此 对 具有 收获 的 
捕食 系统 的 研究 ， 对 生物 资源 进行 合理 的 开发 与 利用 问题 具有 指导 作用 ， 在 经 济 学 和 生物 学 领 
域 都 有 重要 意义 。 目 前 ， 对 于 此 类 问题 的 研究 已 受到 广泛 关注 ， 并 已 取得 一 些 好 的 成 果 [" 引 。 
但 对 于 含有 时 滞 的 捕食 系统 的 研究 ， 很 少见 到 这 方面 的 工作 。 在 捕食 系统 中 ， 捕 食 者 捕食 食 饵 
后 不 是 立即 增加 捕食 者 的 数量 ， 还 需 考虑 捕食 者 应 有 的 妊娠 期 ， 即 捕食 者 捕食 食 饵 后 到 转化 为 
能 量 用 于 繁殖 后 代 所 需要 的 时 间 。 因 此 ， 考 虑 妊娠 期 的 时 滞 因 素 对 系统 的 影响 是 必要 的 ， 更 切 
合 实际 意义 。 时 兆 可 以 将 种 群 稳定 的 正平 衡 位 置 变 成 不 稳定 ， 或 出 现 稳 定性 的 开关 现象 ， 而 有 
些 情 况 下 ， 由 于 时 沿 的 作用 ， 系 统 可 能 会 出 现 非常 数 的 正 周期 解 和 Hopf 分 支 人 6。 

在 文献 了] F, Martin 和 Ruan 在 Wangersky-Cunningham 系统 中 加 入 了 收获 率 ， 得 到 具有 
时 滞 和 常数 收获 率 的 捕食 -被 捕食 系统 


i(t) = z(t)(r1 — az(t) — by(t)) - H, 
Y(t) = —raw(t) + ex(t — 7)y(t — 7), 


其 中 bry 表示 单位 时 间 内 被 y 个 捕食 者 所 吃 掉 的 食 饵 数量 。 从 而 bz 表示 单位 时 间 内 每 一 个 捕食 
者 所 吃 掉 的 食 饵 数量 。 它 除了 与 zx 有 关外 还 反映 了 捕食 者 的 捕食 能 力 ， 称 为 捕食 者 对 食 饵 的 功 
能 性 反应 。 在 系统 (1) 中 ， 功 能 性 反应 被 确定 与 食 饵 的 数量 成 正比 。 比 例 系数 5 标志 着 捕食 能 
力 。 食 饵 数 量 越 大 ， 被 每 个 捕食 者 在 单位 时 间 内 吃 掉 的 就 越 多 。 但 是 ， 捕 食 者 总 有 吃 饱 的 时 
候 ， 因 此 这 种 假定 与 实际 情况 并 不 完全 吻合 ， 功 能 性 反应 与 食 饵 种 群 规模 zx 的 关系 应 该 是 一 种 
包含 饱和 因素 的 函数 关系 。 
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本 文 将 考虑 具有 时 滞 和 Holing- 类 功能 性 反应 的 两 种 群 捕 食 系 统 


aiy(t) ) _ HH, 


i(t) = z(t) (n Sas) — T E ma(f) 


(2) 
azz(t — r)y(t — 7) 
1-4 mz(t-— T) 
其 中 z(t) f y(t) 分 别 表示 食 饵 种 群 和 捕食 者 种 群 在 时 刻 t 时 的 密度 。 参 数 a, ri, r2, al, m, Hi, 
T 均 为 正常 数 ，7 表 示 捕 食 者 种 群 从 幼年 到 成 年 的 成 熟 期 ， 而 且 捕食 者 种 群 只 有 到 成 年 后 
才 具 有 捕食 能 力 ，”i 为 食 饵 的 内 豪 增长 率 ，72 为 捕食 者 的 死亡 率 ，a 为 食 饵 的 密度 制约 系 

35 Hi 为 对 食 饵 种 群 的 常数 收获 率 。 
本 文 我 们 讨论 系统 (2) 将 基于 以 下 初始 条 件 


7(0) = $1(0) 20, y(0)—02(0 20, 8€[-7,0, $1(0)>0, $2(0) >0, 


y(t) 一 r2v(t), 


其 中 
(61(6),02(6)) € C([-7,0], R39), Rio ((9):220,y 20]. 
本 文 的 目的 是 以 滞 量 为 参数 ， 利 用 Hassard 等 在 文献 g] 中 所 介绍 的 规范 型 方法 ， 讨 论 系 
统 (2) 在 正平 衡 点 处 Hopf 分 支 的 存在 性 以 及 分 支 方向 、 分 支 周 期 解 的 稳定 性 ， 最 后 通过 数值 模 
拟 验证 所 得 理论 结果 。 


2 ”正平 衡 点 的 稳定 性 及 Hopf 分 支 的 存在 性 


在 本 节 中 ， 我 们 将 讨论 系统 (2) 正平 衡 点 的 稳定 性 和 Hopf 分 支 的 存在 性 。 
4 z-—ar/n,g —a/ri, £— rat ， 则 系统 (2) 就 变 为 


aD = et-a- A — p 
Pzt — rt = 7) 


l4 yz(t— T) ay(t), 


y(t) = 
其 中 


Ga -—ra/r, f-—as/a, Y= mna, h — aHy/r2, 
为 简便 起 见 ， 仍 以 z, y, tid z, g, to 
Wi —o(8—oy —o)/(B—oy?. RAAH, 4h» hBB—oy—2o» 0 时 ， 系 统 (3) 存 
在 唯一 的 正平 衡 点 E* = (z*, y*)， 其 中 
.a ,AP-a-a) h 


=B ay ” (8 — ay}? a` 


经 计算 可 得 系统 (3) 在 正平 衡 点 E* 的 特征 方程 为 


A? + (pi + p2)À + (—p1À + p3)e *" + pipz = 0, (4) 


其 中 


oy ,aBtoy) h(&—-oy? ao 一 2a 
B BB-a) ag ^ 5 B- oy 


P1220, pa = 
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当 7 = 0 时 ,方程 (4) 即 为 
A? + p2À + pipa 十 Da = 0. 


3492» 0 时 ， 则 由 Hurwitz 判 据 可 知 方程 (5) 的 两 个 根 均 具 有 负 实 部 。 
当 7 > 0 时 ， 注 意 到 iw(w > 0) 是 (4) 的 一 个 根 当 且 仅 当 w 满 足 


w? — p1p2 = —piwsinwr + ps coswT, (pı + p2)w = piw coswT + pa sin wT. 


对 (6) 两 个 方程 分 别 平方 再 相 加 可 得 
wt + paw? + pipi — p — 0. 


H pipo < pa， 方程 (7) 仅 有 一 个 正 根 wo。， 其 中 
w= [3(- pi oi - (tbi - 5) )]-. 


一 Piwl + (pipa + paps * pipa) à 2m dedi d dus 
piug &-- pi wo 
EE, “r= ry (7 = 0,1,…) 时 ， 入 = 土 iwo 是 方程 (4) I] —8 3c 8u £t HEB . 
由 方程 (4) 两 边 同时 对 7 求 导 可 得 


1 ; 
Tj = — arcsin ( 
wo 


dà _ 
[2A + (pi + p2) - (pı T(-mA + p3))e "] i = Xe (—p1 À + p3). 
经 计算 ， 我 们 有 
dA\-! - C ea pj 
Re( 于) l iwo pew +p PS 
从 而 


sf 人 外) eE a) 


根据 文献 [8] 的 第 11 章 的 定理 1.1， 可 得 以 下 结论 。 


定理 2.1 4 wp fi; (j -0,1,2,---) 223! gi (8) 和 (9) 定义 。 如 果 p >0, 2D1D2 < pa» 


() 当 relom) 时 ， 系 统 (3) 的 正平 衡 点 E 是 渐 近 稳定 的 ; 
(i) 当 7 > zo 时 ， 系 统 (3) 的 正平 衡 点 E* 不 稳定 ; 


(üi) Hr =r (7 = 0,1,2,…) 时 ， 系 统 (3) 在 正平 衡 点 B* 附近 出 现 Hopf 分 支 。 


3 ”分 支 方向 及 分 支 周 期 解 的 存在 性 
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(5) 


(6) 


(10) 


(11) 


则 有 


在 本 节 中 ， 我 们 将 用 Hassard 在 文献 上 9 中 所 介绍 的 规范 型 方法 及 中 心 流 形 定理 给 出 系 


统 (3) 的 Hopf 分 支 方向 ， 并 讨论 分 支 周 期 解 的 稳定 性 。 


令 t= sr, z(sr) = zi(s), y(sr) 2 (s), T= To +p, HER, 站 由 (9) 定 义 ， 且 仍 记 上 = s, 


则 (3) 等 价 于 如 下 系统 
| hi) = (m i) (eri - Gl) — 2:0? — Fen)), 


di(t) = (To + u) (Beizi(& — 1) 4- oni (t — 1) — og (t) + Bf (ai (t — 1) gi (t — 1))), 


(12) 
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其 中 
jte Sb Evo) 
97 pae oe ys)" 
cl y'/(0 - vyz*y, c—-1-—2z*—o0, 
did = ($n, 9a)” £ C[-1,0] — C ([-1,0], R?) 定义 一 个 算 子 
L$ = Bé(0) + Có(-1), 


其 中 


sm ( : B eo Je 
0 -a B 


Cl —Q 


且 定 义 


-¢2(0) — f (41 (0), $2(0)) 
BE), p) 三 


f(u, $) = (ro + y) ( 
其 中 了 (81(0), $2(0)) 如 (13) 所 定义 。 


| 
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(13) 


(14) 


于 是 ， 由 Riesz 表示 定理 ， 存 在 分 量 为 上 有 界 变 差 函 数 的 二 阶 和 矩阵 7(9, 4) : [21,0] — R”, 


使 得 对 任意 的 pe C[-1,0|， 有 
0 
i Í GO)dn(0, i). 


事实 上 ， 只 需 取 n(0,y) = Bó(0) — Có(0 -- 1) BlnJ, HH óX Dirac-delta 函数 。 


对 ge Cl[-10， 定 义 
do(0) 
A(u)ó = ó 
ji ansels), — 8-0. 


0 -1<0<0, 
R(u)ó = 
f(n, à), 0 — 0, 


于 是 可 把 方程 (14) 写成 如 下 形式 
Ù; = A(u)U; + R(Q)U, U = (x,y). 
X} y € C1[0,1]， 定 义 
Nit) 0<s<1 
ds ’ Se 
A*y(s) = : 
f eoo. s=0, 
-1 
XT $ e C([-1,0], C?) 1 v e C([0, 1], (C?)*), E XOXZ TERR 


0 0 
(o =B- f [ st mesa, 


(15) 


(16) 


(17) 


(18) 


(19) 


6 和 


其 中 0(9) = 5(6,0, C? 为 复 平面 ， 那 么 算 子 4 = A0) 5A ERAT. nio vim rie 
及 变换 t = st 可 知 土 i7owo 是 算 子 4(0) 的 特征 值 ， 且 其 它 特 征 值 都 具有 严格 负 实 部 ， 从 
而 土 i70wo 也 是 算 子 A 的 特征 值 。 于 是 经 过 计算 ， 可 得 下 面 的 结论 。 

引 理 3.1 q(0) = (1, p)T e078, dr*(s) = D(1,o)T eivoros 分 别 是 A(0) 和 A* 关于 iwoTo, —iwoTo 
的 特征 向 量 ， 而 且 满 足 (q*(s), q(9)) = 1, (a(s),q(0)) = 0， 其 中 


p=P(c—iwo)/a, o=—(c+iwo)e i“ /(8c,), D= (1+ pa + Beo + aape“), 


以 下 ， 我 们 引用 文献 中 中 的 一 些 概念 ， 来 计算 中 心 流 形 Co 在 人 = 0 时 的 坐标 。 令 Ui 是 方 
程 (17) 当 j = 0 时 的 解 ， 定 义 


z(t) = (q*, U), W(t,0) = U.(0) - 2Rez(t)a(0). (20) 
在 中 心 流 形 Co E. 我们 有 W(t,0) = W (z(t), z(t), 0). mE 


2 
W (z, 2,0) = Wao(0) 7 + Wis ()22 + Woo(8- 4, (21) 


2 与 是 中 心 流 形 Co 在 只 和 人 这 两 个 方向 上 的 局 部 坐标 。 如 果 殉 是 实数 ， 则 克也 是 实数 ， 
我 们 仅 考虑 实数 解 的 情况 。 既 然 六 = 0 时 ， 解 V: E Co， 则 与 文献 由 相 同 ， 有 


se) = iuoToz(t) + q" (0) fo(z, Z). (22) 
再 把 (22) 改写 为 
z(t) = iworoz(t) + g(z, Z), (23) 
其 中 
$52) = m5, quii aao rem H (24) 


由 (20)， 我 们 有 
U,(0) = W(z, z, 0) + zq(0) + za(0). 


注意 到 g(0) = (1, p)Teiwom2， 所 以 可 得 
D(z z 
uit(O) = W (z, 2 0) +z 1 e woTo | Zz 1 eiwoTo. 
U2t (8) wo GN Z, 8) p p 


E x i z2 e z z2z 
g(z, Z) = q' fo(z,2) = 9205 十 91122 十 9o2 了 十 921-7 meg (25) 


因此 


第 4 期 田 晓 红 等 :一 类 具 时 滞 和 收获 的 捕食 模型 的 稳定 性 与 Hopf 分支 689 


比较 (24) 和 (25) 的 系数 且 注 意 到 (21)， 可 得 
g20 = 2roD(Gpe 79), gu =270D(—1+Rep(s8 — 1)), 
go2 = 2roD(—5 + 5bpe "ete — 1), 
g21 = 270D | — 67 bpRe{p}r* e70 一 W/L (0) — 88 flppe T 
+ gge (5w c) + plate Piero — pare tiam) 
+ 8 fle "7 (2ppz* — 2z*Re(p) + Ww (-1)) + 68pW (D (—1)e7 «oro 
+ 55 8pWÁS (- 1e" — 2W(? (0) + ppa" — A (0) — p?z* 
+ 2pz* — WA? (0) + gz* — Wf (0)p — 5 WÍ O) + 6Re(p)]. 
现在 我 们 计算 Wao(0) 和 Wi1(9)。 根 据 (17) 81 (20) 我 们 有 


AW — 2Re(q* (0) foa(0)}, -1€8«0, 


W = U-iq-zd- 
AW — 2Re(q*(0)foq(0)) + fo, 0 — 0, 


def 


kí AW + H(z,z,0), 
其 中 : n 
H(z,2,0) = Ha. 4 Hiizz 十 Ha Rossi (26) 
当 0 e [-1,0) 时 ， 我 们 可 以 得 到 
(A — 2iworo) Wa0(0) = —Hao(6),  AWai (6) = —Hai (0). (27) 
BT 
H(z, 2,0) = —2Re(q'(0)foq(8)) = —9a(0) — ga(), (28) 
比较 (28) 两 边 的 系数 可 得 
Hao(0) = —9209(0) 一 5025(0)， Hi1(0) = —9114(8) — 211a (8). (29) 
另 一 方面 ， 我 们 有 
Wao(8) = 2i woToW20(0) + g209(0) 十 5025(0). (30) 


注意 到 g(9) = (1, p) Tier, Re 


ig20 iwoTo00 igo2 = —iwoToO 2iwoT00 
0 = 0TO 070 woTo T 
Wao(0) 2s dU 十 c5 e + Eje (31) 
类 似 可 证 
Wi1(0) = -1 dg(0)eiwonne + 911 5(0)erieone + pr. (32) 
WoTO WoTo 
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EH Ei, Ez 的 值 可 由 0 = 0 时 五 的 值 来 确定 。 由 A4 的 定义 和 (27) 可 得 
i : dn(0)Wao(0) = 2iworoW20(0) — Hao(0), J dn(6)W11 (0) = -H11 (0), (33) 
XÆ n(0) = $(0,0). WHE H KEX, RITE 
H20(0) = —g204(0) 一 5025(0) + T041, | Hii(0) = —g114(0) — 9119(0) + T042, 
其 中 
Ai = (—2 - 2p, 28pe P9), Az = (—2 — 2Re(p), 26Refoh) . (34) 
我 们 将 (31) 和 (34) 代入 到 (33)， 由 此 解 得 


1 
B= $ (AL, AD, 


A 

其 中 

Al = —4uj — ca + 2(a — cjwoi + (c + c — 2iwo)are o, 

Al = —2o(1- p) — 4wo(1 + p)i + 2œe7 %00, 

A? = 2Be ?vo (2pwoi — cp — (1 + p)c1), 
同 理 可 求 出 

E; = 元 (Ab ABT, 

其 中 


As = acl， Al = —2aRe{p}, A3 = 26Refp}(cel — c) 28e. 


通过 上 面 的 分 析 ， 我 们 可 解 出 Wao(0) 和 Wi1(9)。 更 进一步 , 我 们 可 以 求 出 gz1, 那么 (25) 中 
的 gi; 就 完全 可 以 由 系统 (12) 的 系数 表 出 了 。 于 是 ， 由 文献 [10] 中 的 一 般 性 定理 可 以 计算 出 下 
列 数值 


C1(0) = xoc (auem - 2|gul? 一 5lgoz) 十 i 
wa 人 p, = Ref (0), (35) 
q, — Inte (0)) + palm{X (ro)} 

Te . 


woTo 


Bi (35) 中 C,(0) 的 表达 式 ， 易 得 出 12, bo T2 的 值 。 另 一 方面 ， 我 们 知道 u 决定 了 Hopf Zt xx 
的 方向 ， 如 果 ji > 0(< 0)， 则 Hopf 分 支 是 超 临 界 ( 亚 临 界 ) 的 。 当 r > m (< Tro) 时 ，Hopf 分 
支 的 周期 解 存在 ，0p2 决定 了 分 支 周 期 解 的 稳定 性 : 250, < 0(> 0) 分 支 周期 解 是 稳定 的 (不 稳 
ER) To 决定 了 分 支 周 期 解 的 周期 ， 如 果 T2 > 0(< 0) 周期 是 增加 (减少 ) 的 。 
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4 ”数值 模拟 
本 节选 取 a = 1, 8 =4, 7=1.5 h = 0.1， 则 系统 (3) 变 为 
=a- AO T 
; 4z(t — r)y(t — 7) G6) 
WO = a IO 
计算 可 得 


To =œ 0.0986, C1(0) ~ —0.4212 + 0.2714i. 


当 7 = 时 ，fB。< 0， 周 期 解 是 稳定 的 ;jpa = 5.1024 > 0，Hopf 分 支 是 超 临界 的 ; T, = 
1.7460 > 0， 周 期 是 递增 的 。 当 r < no 时 系统 (36) 的 正平 衡 点 是 渐 近 稳定 的 ， 如 图 1; 
当 7 > To 时 系统 (3) 的 正平 衡 点 是 不 稳定 的 ， 如 图 2。 


0.5 
x ~ 
5 5 «| 
G EA 
$ $ 
ES 5 03 
0.2 
d m 0.1 
0 200 400 600 2 4 6 
t prey x 
图 1: Rr = 0.02 < m， 正 平衡 点 局 部 渐 近 稳定 
0.8 — 一 一 1 1 -一 一 
| | 
06 0.8 | | 0.8 
x 5 06 | 5 06 
全 
504 5 0.4 
0.2 
| 0.2 0.2 
0 0 à 0 UN E 
0 100 200 300 0 100 200 300 0 0.2 0.4 0.6 0.8 
t t prey x 


图 2: 取 T 二 0.15 > mm， 正平 衡 点 不 稳定 
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Stability and Hopf Bifurcation of a Predator-prey Model with Time 
Delay and Harvesting 


TIAN Xiao-hong, XU Rui, WANG Lili 


(Institute of Applied Mathematics, Mechanical Engineering College, Shijiazhuang 050003) 


Abstract: A predator-prey model with harvesting and a time delay describing the gestation period of 
the predator is considered. We investigate the local stability of a positive equilibrium by analyzing the 
corresponding characteristic equation. It is proved that when the time delay 7 increases, the positive 
equilibrium looses its stability and a Hopf bifurcation occurs when 7 passes through a critical value. 
Formula are derived to determine the direction of bifurcations and the stability of bifurcating periodic 
solutions by using the normal form theory and the center manifold theorem. Numerical simulations 
are carried out to illustrate theoretical results. 

Keywords: predator-prey system; time delay; functional response; Hopf bifurcation; stability 
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